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Abstract. The paper is devoted to the range description of the 
Radon type transform that averages a function over all spheres 
centered on a given sphere. Such transforms arise naturally in 
thermoacoustic tomography, a novel method of medical imaging. 
Range descriptions have recently been obtained for such trans- 
forms, and consisted of smoothness and support conditions, mo- 
ment conditions, and some additional orthogonality conditions of 
spectral nature. It has been noticed that in odd dimensions, sur- 
prisingly, the moment conditions are superfluous and can be elim- 
inated. It is shown in this text that in fact the same happens in 
any dimension. 



This papers deals with the integral transform that associates to a 
function / G C(IR n ) its spherical means 



where dO is the normalized area measure on the unit sphere 



This integral geometric transform arises in a variety of pure and ap- 
plied mathematics topics from PDEs, medical imaging, and geophysics 
[3, 4, 8, 9, 11, 18, 23, 24], which we will not survey here. 

Experience of integral geometry and its applications (e.g., [10, 14- 
17,27]) shows importance of the range descriptions for integral geo- 
metric Radon type operators. For the spherical mean transform / i— > g 
introduced above, the range is provided by the famous Asgeirsson's 
theorem [7, 9, 17, 18], which says that the resulting functions g(x, t) are 
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exactly the solutions of the following initial value problem in IR n for 
the Darboux-Poisson-Euler equation: 

(i) <g(x,o) = f(x) 

[§f(*,o) = o. 

Note, however, that the spherical mean operator / i— > g increases 
the number of variables from n to n + 1. It is thus natural to reduce 
the choice of spheres to an n-dimensional family. This can be done 
in various ways. One is to allow only a discrete set of values of the 
radius t. This kind of a restricted spherical mean transforms has been 
intensively studied for many years in harmonic analysis (see [32, 36] 
and references therein) and will not be discussed here. Our target 
transform is different, the one where the set of centers x is restricted 
to an {n — l)-dimensional variety S C IR n . This restricted type of the 
operator was studied, due to its relevance to approximation theory, in 
[4,23,24]. We in fact restrict this transform even further, assuming S 
to be the unit sphere 

S:={xe R n : \x\ = 1}. 

Moreover, it will be assumed that the original function / is supported 
in the unit ball 

B := {x E W 1 : \x\ < 1}. 

Thus, our article is devoted to the study of the following restricted 
spherical mean operator Rs'. 

g(x,t) := R s f(x,t) := f(x + t9)d6, (x, t) G S x R+. 

J\e\=i 

The transform just introduced has been studied intensively in recent 
years (see surveys [2,13,20,31]), since it represents the data collected 
in, and thus is the major component of the newly developing method 
of medical imaging called thermoacoustic tomography [34,35]. 

Our interest here is in obtaining the range descriptions (i.e., analogs 
of Paley- Wiener theorem) for the operator Rs on the space C£°(B). 

It is instructional to look at the example of the standard Radon 
transform R that produces integrals g(9, s) of f(x) over (n— l)-dimensional 
planes xx6 = s. The well known [10, 14, 17, 27] range description for R 
contains (besides some natural smoothness requirements) three types 
of conditions: 

• the natural support condition on g = Rf: planes not touching 
the support of / produce zero integrals, 
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• the so called moment conditions: Mk(6) := f s k g{6,s)ds is the 
restriction to the unit sphere of a homogeneous polynomial of 
degree k, and 

• the evenness condition g(6, s) = g(—0, —s) that reflect the fact 
that the planes given by the equations x-6 = s and x-(—9) = —s 
coincide. 

On the face of it, the moment conditions constitute the only non-trivial 
part. However, this is where the standard Radon transform is mislead- 
ing. Indeed, for various its generalizations the situation changes. For 
instance, for some weighted (so called attenuated) Radon transforms 
arising in emission tomography, the analogs of moment conditions were 
found rather fast, while it took many years for an analog of the "triv- 
ial" evenness condition to be found [21,28]. A similar thing has hap- 
pened with the transform R$: moment type conditions were explic- 
itly described in [30], although they were implicitly present already in 
[4,23,24]. It was harder to find the complete set of conditions, i.e. 
to find analogs of the evenness condition. This was done in 2d in [5], 
then in odd dimensions in [12], and finally various range descriptions 
that apply to any dimension were obtained in [3] . In most cases, range 
descriptions were in a sense parallel to the Radon case: they consisted 
of a natural support condition for g = Rsf, the moment conditions, 
and some "orthogonality" condition that thus looked like an analog of 
evenness in the Radon case. However, the following peculiar thing was 
observed in odd dimensions in [12]: the moment conditions could be 
dropped completely from the range description. In other words, the 
orthogonality conditions were apparently not just analogs of evenness, 
but also contained the moment conditions as well. This absolutely 
non-transparent fact was later confirmed in the odd dimensions for the 
different range descriptions obtained in [3]. The initial impression of 
the authors was that the situation might depend on the validity of 
the Huygens' principle for the wave equation, which only holds in odd 
dimensions [9]. However, the goal of this text is to show that the mo- 
ment conditions can be dropped from the range descriptions 
in any, not only odd, dimension. 

The next Section 2 contains the formulation of the result, which 
is then proven in Section 3. The paper ends with the Remarks and 
Acknowledgments sections and the bibliography. 

2. Formulation of the result 
Let us first introduce some notations. 
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We use the standard notation J P (X) for Bessel functions of the first 
kind. We will also need to use what is sometimes called normalized or 
spherical Bessel functions [26] : 

2*T(p+l)J p (A) 



(2) J P (A) = 



For p = according to Poisson's representations of Bessel functions 
[18, 26], this is the spherical average of a plane wave in IR n . 

Consider for each m G Z + an orthonormal basis Y™, 1 < I < d(m) 
in the space of all spherical harmonics of degree m in 1". 

Let g(x, t) be a function on S x R + . This function can be expanded 
into the Fourier series with respect to spherical harmonics as follows: 

(3) 0(s > O=X>m > l(*)*r(s), 

m,l 

where 

(4) g nhl (t)= [ g(6,t)YT(e)d0. 

Jees 

We will also need the Fourier-Bessel transform of g m j. 

POO 

9m,l( X ) = / gmMinn-l^t^dt. 

Jo 

The main result of this paper is the following: 

Theorem 1. The following three statements for a function g G C^(Sx 
[0,2]) are equivalent: 

(1) The function g is representable as Rsf for some f G C™(B). 

(2) For any eigenvalue —A 2 of the Dirichlet Laplacian in B and the 
corresponding eigenfunction i])\{x), the following orthogonality 
condition is satisfied: 

(5) J g(x,t)d^ x (x) Jn/2 _ 1 (\t)t n - 1 dxdt = 0. 

Sx[0,2] 

Here d u is the exterior normal derivative at the boundary of B. 

(3) For any integer m and 1 < I < d(m), function g mt i(X) van- 
ishes at all non-zero zeros of the Bessel function J m+n / 2 -i(A) 
(equivalently, at all zeros of j m + n /2-i(X)). 

Not all the claims of this Theorem are new. Indeed, implication 
(1) =3- (2), (3), as well as equivalence (2) <^ (3) were established in 
[3, Theorem 10]. In fact, it was also shown there that each of the 
statements (2) and (3) implies (1), if being supplemented by additional 
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Moment conditions:: for any non-negative integer k, the func- 
tion M k on S defined as 
2 

(6) M k {x) := J t lk+n - 1 g{x,t)dt, x G S 

o 

has an extension to W 1 as a polynomial of degree at most 2k. 
Therefore, due to [3, Theorem 10], the proof of Theorem 1 will be 
concluded, if we prove the following surprising result: 

Theorem 2. For a function g G C^°(S x [0, 2}), any of the equivalent 
conditions (2) and (3) implies the moment conditions (6). 

3. Proof of Theorem 2 

Since the orthogonality condition 3 of Theorem 1 is formulated in 
the Fourier domain, we start the proof with reformulating the moment 
conditions in a similar way, as it was done in [3]. The next lemma is 
proven in [3, Lemma 9]. For completeness, we provide its brief proof 
here. 

Lemma 3. [3, Lemma 9]) The following conditions on a function g G 
Cq 3 (S x [0,2]) are equivalent: 

(1) The moment conditions (6) for g(x,t). 

(2) For any m > 1 and 1 < I < d{m), function g nh i{\) vanishes at 
A = to the order at least m. 

Proof. Expanding g(9, t) into an the orthonormal basis Y™ of spherical 
harmonics on S, one obtains 

oo oo 

(7) M k (9) = J t^-'giO^dt = J2 Y i m ( ) J t^-'gmMdt. 



A spherical harmonic Y™ of degree m can be extended to a polynomial 
of degree d if and only if d > m and d — m is even. Thus, the moment 
conditions require that 



I 



oo 

t 2k+n - 1 g m ,i(t)dt = 

'o 

for m > 2k. 

On the other hand, using the well known series expansion j p (X) = 
YlT=o CfcA 2fe for Bessel functions, one arrives to the formula 



oo 



9m,l(X) 



J2 C k^ 2k / t 2k+n - x I Y l m (6)g(6,t)dS(e)dt. 



k ~° o ees 
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This shows that the moment conditions are equivalent to the vanishing 
of all terms in this series with 2k < m. □ 

We now return to the proof of Theorem 2. We assume that g G 
C^°(S x [0, 2]) satisfies conditions 3 of Theorem 1, i.e. that for any m > 
1, 1 < Z < d(m), function <? m ,;(A) vanishes at all zeros of j m + n /2-i(ty ■ 
We need to show that then A = is a zero of order at least m oig m j(\). 

In fact, we will achieve this working with each individual function 
{Jm,i(ty- To simplify notations, from now on, when considering a func- 
tions g m j or g m j, we will drop in the notations the irrelevant index /, 
which does not enter any considerations. Thus, g m will denote any of 
the functions g mj i with 1 < I < d{m). 

So, let g satisfy the conditions above. Then g m (A) is an entire func- 
tion of the following Paley- Wiener class: 

(9) |? m (A)| < C N (1 + |A|)-VI 7mA ' for any N > 0. 

We also know that it vanishes at all zeroes of the entire function 
im+n/2-i(A), which happen to be simple. Thus, the ratio H(\) of these 
two functions is also entire: 

(10) g m (X)=Jm + n/2-lWH(X). 

We will now determine the Paley- Wiener class to which H belongs. 
This is done using the following simple (and probably well known) 
estimate from below for the Bessel function of the first kind J v : 

Lemma 4. [5] On the entire complex plane, except for a disk Sq cen- 
tered at the origin and a countable number of disks Sk of radii n/6 
centered at points n(k + ^j^), one has 

(11) |j w (z)| > c>0 - 

For completeness, we present here the proof of this lemma, borrowed 
from [5]. 

Proof of Lemma 4: We split the complex plane into three parts 
by a circle Sq of a radius R (to be chosen later) centered at the origin 
and a planar strip {z = x + iy\ \y\ < a}. Here part I consists of points 
z such that \z\ > R and \Imz\ > a; in part II, \z\ > R and \Imz\ < a; 
finally, part III is the interior of So, i.e. contains points z such that 
\z\ < R. It is sufficient to prove the estimate (11) in the first two parts: 
outside and inside the strip. In fact, it suffices to consider only the 
right half plane Re z > intersected with parts I and II. 

The Bessel function J v {z) has the following known asymptotic rep- 
resentation in the sector |arg^| < n — 5 (e.g., [6, formula (4.8.5)] or 
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[22, formula (5.11.6)]): 

Ju{z) = a/^cos^ - f - f )(1 + 0(\zn) 



(12) 

-y^^-f -f) (T + o(kr 3 )) 

Let us estimate it in the part I, where \Imz\ > a and \z\ > R for 
sufficiently large a and R (and, as we have agreed, Rez > 0). There 
sin 

one has = cos z (Od^ -1 )), and thus (12) implies 



For sufficiently large a, i? this leads to the needed estimate 

/7 p |Jm*| 

(13) \Uz)\> ' 

In the part 11 of the plane (right half of the strip), due to boundedness 
of sinz, one has 

Mz) = ^ [cos(, - ^ - J)(i + o(N~ 2 )) + od^r 1 ) 

Let us consider the system of non-intersecting circles Sk of radii | 
centered at Zfc = | + kn + ^ + |. Outside these circles, | cos(,2 — ^ — 
f)|>C and thus 

i^i^a+od^r 1 )), 

which implies, inside of the strip and outside the circles Sk, the estimate 

n e \imz\ 

(14) \J u (z)\ > 

for | z | > R with a suitably chosen sufficiently large R. This proves the 
statement of the lemma. □ 

Corollary 5. Function H(\) in (10) is even, entire, and satisfies the 
following Paley- Wiener estimate 

\H(X)\ < C N {1 + |A|) _Jv e |/mA| /or any N > 0. 

Proof. Since H = g m / jm+n/2-1, this follows from (9) and (11). □ 

The well known Paley- Wiener theorem for the Fourier-Bessel trans- 
form (see, e.g. [3, Lemma 1] and references in that paper) shows that 
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then the following representation is possible: 

POO 

(15) H(X) = / h(t) 3(n _ 2)/2 (t)t n ~ l dt, 

Jo 

where h(t) G Cq°(R) is an even function on the real line supported in 
the interval [—1,1]. 

We introduce now the function h on IR n as 



%) = H\y\), where \y\ = y J^J/J. 

Then h is clearly a C^-function on M. n supported in the unit ball -8(0, 1) 
centered at the origin. 

We would like to lift similarly the Bessel function j m + n /2-i to W 1 . 
This produces the function 

E m (v) := j m +n/2-l(\v\) 

Lemma 6. Function E m {y) is, up to a constant factor, the n- dimensional 
Fourier transform of the function (1 — \y\ 2 )'"^~ 1 ,y G W 1 , i.e. 



s 



(16) E m {v) = const / e- w »(l - \y\X dy 
Here u + := max(u, 0). 

Proof. This identity follows from the Sonine integral formula [6, for- 
mula 4.11.11], but we will indicate the proof, to make the presentation 
more self-contained. 

The integral on the right hand side is clearly invariant under ro- 
tations, so we may suppose that v = \v\ei. We set y = (s,y') with 
y' G R n_1 . Applying Fubini's theorem, the integral on the right can 
then be written 

1 e-' lMs [ (1-s 2 - \y'\ 2 )^- 1 dy' ds. 

Js=\v\ 

The substitution y' = y/l — s 2 a reduces this to 

withC m> „ = l\S n - 2 \B(v=±,f), w here \S n ~ 2 \ is the area of the (n-2)- 
sphere, and B(-,-) is Euler's Beta function. The integral is a constant 
multiple of a well known representation of j m + n /2-i- 1=1 

Corollary 7. Function g m (\v |) on M. n is representable as the convolu- 
tion 

(17) g m (\v\) =w m (\v\)*h, 
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where w m { \v\) = (1 — \v J 2 )™" 1 and h G C^°(IR n ) is supported in the unit 
ball B {0,1). 

Indeed, this follows from (10) and Lemma 6. 

Let us now pay attention to the behavior of the function g m at the 
point t — 0. Namely, since g m (t) G C°°[0, 2], it vanishes at t — with 
all its derivatives. This behavior is inherited by the smooth function 
<7 m (|i;|) in W 1 : it vanishes at the origin with all its derivatives. In 
particular, if D — (d/dvi, ...d/dv n ) and a = (cci, a n ) is a multi- 
index, we get the identity 



(18) = D a g m (0) = w m * (D°h) \ v=0 = J (l-M 2 )?" 1 (^Hv)) 



dy. 



Since h is smooth and supported in the unit ball, the right hand side 
is equal to 



J (I ~ M 2 )™- 1 (D a h(y)) dy. 



Integrating by parts, we now conclude that the following equalities hold 
for any a: 

(19) J (D°(l-\y\ 2 ) m - 1 )h(y)dy = 0. 

We can now establish the claim central to the proof of Theorem 2: 
Lemma 8. For any k < m — 1, the following moment relation holds: 

(20) J \y\ 2h h(y)dy = J \y\ 2h h(y)dy = 0. 

1" \y\<l 

Proof. Indeed, expanding (1 — |y| 2 ) m_1 according to the binomial for- 
mula and using (19), leads in particular to 

E (!)(-!)* / ^(\y\ 2k )Hy)dy = o, 

!1 |y|<i 

for any integer j > 0, where A J is the j-th power of the n-dimensional 
Laplace operator. Using the well known and easily verifiable relation 
A\y\ 2k = ct\y\ 2k ~ 2 with Ck ^ 0, and running the value j from j = to 
j — m — 1, one obtains a uniquely solvable triangular homogeneous lin- 
ear system of equations with respect to the moments f (\y\ 2k )h(y)dy 
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for < k < m — 1. Thus, these moments are all equal to zero. Due to 
the evenness of h, all odd order moments also vanish. □ 

The vanishing of the moments shown in the last Lemma implies 
that the Fourier transform h of h, and hence the corresponding one- 
dimensional function H, has zero at the origin of order at least m. 
Due to (10), this feature is inherited by g m: which finishes the proof of 
Theorem 6, and thus also of Theorem 1. 

4. Remarks 

The implication (1) in Theorem 1 =>- moment conditions was 
implicitly present in [4, 23, 24] and explicitly formulated as range 
conditions in [30]. 

If S is not a sphere, most probably the moment conditions are 
independent of the orthogonality ones. 

We have shown that in fact g m) i has at the origin a zero of the 
order 2m, not just mere m. As it was mentioned in a similar 
situation in [3], this feature reflects the fact that S is a sphere 
and most probably will not hold for more general surfaces S. 
The range conditions provided in [3] had also another incarna- 
tion, besides the versions (2) and (3) of Theorem 1. Indeed, 
it was shown in [3, Theorem 3] that the following condition is 
equivalent to (2) and (3): 

The backword in time solution of the Darboux equation (1) with 
the boundary data provided by the function g and zero initial 
conditions for a sufficiently large moment T of time does not 
develop a singularity when t — > 0. (see details in [3]). 

This condition can be substituted instead of (2) or (3) in 
Theorem 1. 
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